Link-disorder fluctuation effects on synchronization in random networks 
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We consider one typical system of oscillators coupled through disordered link configurations in 
networks, i.e., a finite population of coupled phase oscillators with distributed intrinsic frequencies on 
a random network. We investigate collective synchronization behavior, paying particular attention 
to link-disorder fluctuation effects on the synchronization transition and its finite-size scaling (FSS). 
Extensive numerical simulations as well as the mean-field analysis have been performed. We find that 
link-disorder fluctuations effectively induce uncorrelated random fluctuations in frequency, resulting 
in the FSS exponent v — 5/2, which is identical to that in the globally coupled case (no link disorder) 
with frequency-disorder fluctuations. 
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I. INTRODUCTION 

Recently, it has been reported that disorder fluctua- 
tions play a crucial role on the critical scaling near the 
onset of synchronization, even in the mean-field (MF) 
regime [lr|3|. In particular, one recent study on the col- 
lective synchronization in the system of globally coupled 
oscillators (Kuramoto model) reports that "sample-to- 
sample fluctuations" induced by different realizations of 
intrinsic random frequencies yield significantly different 
finite-size scaling (FSS) from the ordinary MF predic- 
tion |l|. The anomalous fluctuations are characterized 
by the FSS exponent v — 5/2 [l|. When the random dis- 
order in distributed intrinsic frequencies is completely re- 
moved, fluctuations are still anomalous but much weaker 
(p ~ 5/4) [2J, because dynamic (temporal) fluctuations 
dominate in this case. We call the latter as the noise- 
less frequency distribution, while the former as the noisy 
frequency distribution. 

When the system is put on a complex network, we 
need to consider another kind of fluctuations, i.e., "link- 
disorder fluctuations" induced by different configurations 
of connectivity among vertices in the network. Consider- 
ing that many systems in nature are based on a complex 
network topology, it would be interesting to investigate 
how the link-disorder fluctuations affect the collective be- 
havior of a given system in general. Most of collective 
behaviors invoking phase transitions in networks belong 
to the conventional MF class where dynamic (thermal) 
MF fluctuations are dominant, thus sample-to-sample 
disorder fluctuations become usually irrelevant |l|, 14-11 lj . 
However, in the presence of additional distributed quanti- 
ties at vertices, link-disorder fluctuations would generate 
effective disorder fluctuations in distributed quantities, 
which might yield non-MF fluctuations. 

In this study, we consider the system of coupled oscil- 
lators on a sparse link-disordered network. To separate 
out the effects of uncorrelated random link-disorder fluc- 



tuations, we focus on the system with a noiseless intrinsic 
frequency distribution on the Erdos-Rcnyi (ER) randonr 
and z-regular network [12j. The FSS property is inves- 
tigated through the phase synchronization order param- 
eter in extensive numerical simulations and the value of 
the FSS exponent v is estimated. We find that v ~ 5/2, 
which is identical to that of the globally coupled model 
with the noisy frequency distribution. This implies that 
the random link-disorder fluctuations effectively generate 
random and independent disorder in the frequency distri- 
bution. It is also found that the addition of random dis- 
order in the intrinsic frequency distribution (noisy case) 
on the ER network does not change the FSS exponent. 

The present paper is outlined as follows. In Sec. II we 
introduce the model, and explain how we generate the 
noiseless frequency distribution. In Sec. Ill, we present 
the mean-field (MF) analysis on the ER network with the 
predictions on the order parameter exponent j3 and the 
FSS exponent v. Section IV shows numerical simulation 
results. A brief summary is given in Sec. V. 



II. MODEL 

We begin with a finite population of N coupled phase 
oscillators on a sparse network. To each vertex i of the 
network, we associate an oscillator whose state is de- 
scribed by the phase <f>i governed by the equation of mo- 
tion 
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= Ui-K^ 



an sm 



~<t>j), 



(1) 



where Wj represents the intrinsic frequency of the ith os- 
cillator. It is assumed that {tOi} are distributed according 
to a symmetric and unimodal distribution function g(co) 
such as a Gaussian one with zero mean ((w) = 0) and 
finite variance ((ujiuij) = a 2 SiSj). The {a^-} denotes the 



adjacency matrix denned by 



1, when i and j are linked, 
0, otherwise, 



(2) 



and the degree ki of the vertex i is defined as the num- 
ber of vertices linked to i, i.e. ki = J2j a ij- I n a sparse 
network, the total number of links is proportional to N. 
The second term on the right-hand side of Eq. (p} rep- 
resents ferromagnetic coupling to neighboring oscillators 
on the network (K > 0), so the neighboring oscillators 
favor their phase difference minimized. 

When the coupling is weak (small K), each oscilla- 
tor tends to evolve with its own dynamics described by 
dcj)i/dt ~ u)i, showing no synchrony. As the coupling 
is increased, the oscillators interact with each other and 
eventually a macroscopic number of synchronized oscilla- 
tors can appear if the coupling is strong enough. Emer- 
gence of macroscopic synchronization can be probed by 
the phase order parameter [13( 
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where A measures the magnitude of the phase synchrony 
and 9 denotes the average phase. 

Near the synchronization transition (A « 0), the MF 
analysis shows that the order parameter A vanishes with 
the ordinary MF exponent (3 = 1/2, if the degree distri- 
bution {ki} is not strongly heterogeneous [5|, |6|. The ER 
random network has an exponential degree distribution 
with weak heterogeneity, thus the ordinary MF exponent 
should be found. 

The FSS is governed by fluctuations. In the system 
described by Eq. ([1]), there could be three different fluc- 
tuations. One comes from dynamic (temporal) fluctua- 
tions which persist even in the steady state due to the 
presence of running (non-static) oscillators. Two oth- 
ers are due to quenched disorder in link configurations 
{a^} in a sparse network and possibly in frequency dis- 
tributions {uJi} obtained by a random drawing procedure 
from the distribution g(uj). It is difficult to study ana- 
lytically the dynamic fluctuations even in the globally 
coupled case (no link disorder) and there still exists a 
controversy [2|, |3[. However, it is known that frequency- 
disorder fluctuations dominate over dynamic fluctuations 
at least in the globally coupled case [l|, [2| . 

In this work, we focus on the effects of link-disorder 
fluctuations separately, so that we need to remove 
frequency-disorder fluctuations completely from the sys- 
tem. This can be achieved by generating the frequencies 
{u>i} following a deterministic procedure given by [2J 



i-0.5 

N 



g{uj)du). 



(4) 



The generated frequencies are quasi- uniformly spaced in 
accordance with the distribution g(co) in the population. 



As this noiseless frequency set is uniquely determined, 
there is no disorder originating from different realiza- 
tions of frequencies. For the globally coupled model with 
this noiseless set, dynamic fluctuations dominate the FSS 
with the weak FSS exponent v ~ 5/4 [2j , in contrast that 
v = 5/2 [l| with disorder in frequency sets by randomly 
and independently drawing frequencies from the distri- 
bution g(u>). 

The question we raise in this work is what kind of FSS 
emerges in a sparse network such as the ER network with 
the noiseless frequency set, i.e. what is the role of disorder 
in link configurations? Also does the FSS change or not 
if disorder in frequency sets is added in the ER network? 
In the following sections, we start with a MF analysis 
and report extensive numerical results. 



III. MEAN-FIELD ANALYSIS 

We consider the ER networks as simple examples of 
uncorrelated random sparse networks. In the ER random 
network, any pair of vertices is linked independently with 
probability z/N. In the large N limit, the average degree 
becomes (k) = z and the de gree distribution is given by 
the Poisson distribution [12|, |l4 1 



P(k) 
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(5) 



In the z-regular random network, all vertices have the 
same degree z with random connectivity between ver- 
tices, so the degree distribution is given by the Kronecker 
(5-function as P(k) = Sk, z - 

Taking the annealed MF approximation for connec- 
tivity of the network a^ « kikj/N(k) (heterogeneous 
MF theory), one can introduce the global ordering field 
Hexp{i6) Q as 
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and then the equation of motion, Eq. (TTJ) , can be rewrit- 
ten as 



u>. 



k;KH sinf 



(7) 



In the long-time limit where the global ordering field 
approaches a constant in average, one can easily solve 
Eq. ([7} for each fa with a constant H and 9. Then, we 
can set up the self-consistency equation for H through 
Eq. © as 



A 



H N^(k)y l \k 3 KH 
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kjKH 



(8) 



where Q(x) is the Hcavisidc step function, which takes 
the value 1 for x > and otherwise. Each term inside 



the summation represents the contribution from an oscil- 
lator with uij at vertex j with degree kj . Note that only 
entrained oscillators with \uij\ < kjKH contribute. 

The frequency set {iOi} fluctuates randomly over the 
samples for the noisy frequency distribution, and the de- 
gree configuration {ki} also fluctuates randomly for the 
ER random network. Thus, the nontrivial (nonzero H) 
solution of the self-consistent equation also fluctuates 
over the samples. With either one of fluctuations, we 
expect the Gaussian fluctuations for sufficiently large iV 
from the central limit thcorem,and the self-consistency 
equation for small H becomes [6[ 

H = (K/K C )H - c{KHf + d{KH) l/2 N- l/2 (, (9) 

with K c = [2/Trg(0)](k)/(k 2 ), and constants c = 
-[^"(0)/16](fc 4 )/(fc) andd = v /[4 5 (0)/3](fc 3 )/(fc 2 ). The 
term £ is a Gaussian random variable with zero mean 
and unit variance, which represents frequency-disorder 
fluctuations and/or link-disorder fluctuations. 

The scaling solution to Eq. ([9]) can be easily obtained 
in the FSS form [l,[lH 



H = N~ 1/5 f[(K - K C )N 



2/5i 



(10) 



implying f}/v = 1/5 and v = 5/2. We note that this 
result is the same as that of the globally coupled noisy 
oscillators with frequency-disorder fluctuations |f|, l2| ■ 

One exceptional case can be constructed when the sys- 
tem is put on the z-regular random network with the 
noiseless frequency set. In the present MF scheme, there 
is no fluctuation in the solution of H and one may expect 
the same FSS behavior as that of the globally coupled 
noiseless oscillators with /3/p ss 2/5 and v ps 5/4, rep- 
resenting dynamic fluctuations only. However, the MF 
approximations ignore fluctuations in connectivity disor- 
der (absent for the globally coupled case), which might 
give rise to another type of sample-to-sample Gaussian 
fluctuations dominant over dynamic fluctuations. 

In the next section, we check the validity of the MF 
analysis via extensive numerical simulations in the ER 
random and z-regular networks with the noiseless and 
noisy frequency distributions. 



IV. NUMERICAL ANALYSIS 

We perform extensive numerical simulations for the 
system governed by Eq. (JTJ. The ER and the z-regular 
random networks are generated for the average degree 
z = (k) = 6 up to the system size TV = 12800. For 
the noiseless frequency distribution, we generate the 
intrinsic frequencies by the process shown in Eq. (|4"|) 
with g(u>) = (27T(t)~ 1 / 2 exp(— uj 2 /2a 2 ) with unit variance 
{a 2 = 1). We used Heun's method [16[, with a dis- 
crete time step St = 0.01, in the numerical integration 
of Eq. jl} up to 4 x I0 4 time steps. Initial values for 
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FIG. 1: (Color online) (a) The order parameter A is plotted 
as a function of the coupling strength K for various system 
size TV. (b) Binder cumulant T3a versus K for various system 
size TV, where we estimate K c — 0.300(5) from the crossing 
points. 



{4>i(0)} are chosen random, and data are collected only 
for the latter half of time steps to avoid any transient 
behavior [171 ]. For each network size, we also average the 
data over 10 2 ~ 10 3 different realizations of the network. 
First, we consider the ER random network and mea- 
sure the phase order parameter A as a function of the 
coupling constant K for various system size TV, see 
Fig. H] (a). The critical coupling strength K c , beyond 
which the phase synchronization occurs, can be accu- 
rately located by utilizing the Binder cumulant [l8| 



B A = 1- 



A 4 
3A 2 ' 



(11) 



where the overbar represents the time average after suf- 
ficient relaxation and [• • • ] represents the average over 
different realizations of the networks, respectively. In 
the weak coupling regime K <C K Cl the Binder cumulant 
T3a goes to 1/3 in the large TV limit, while it goes to 2/3 
for the strong coupling regime (K 3> K c ). The Binder 
cumulant curves with different TV cross each other at the 
critical coupling strength K c , as seen in Figure Q] (b). 



We estimate K c = 0.300(5) by extrapolating the crossing 
points in the large N limit. Note that the MF prediction 
on K c is about 0.23 in Eq. ©, which is quite far away 
from the numerical estimate. It is not surprising to see a 
considerable shift of the critical point, because the inher- 
ent linking disorder in the quenched networks generates 
finite correlations in neighboring nodes. However, it is 
hoped that the scaling nature near the transition is not 
affected by finite correlations, since is is usually universal. 
We now check the FSS relation of Eq. (fT0| . As one can 
easily show A ex H for small H at the MF level [6j, we 
expect 

AiV 1/5 = f[{K - K c )N 2/ % (12) 

where the scaling function f[x] is characterized by 



f[x] ~ { x 1 ' 2 , x > 0, (13) 
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for the ER random networks with the noiseless frequency 
distribution. Figure [2] shows the scaling function f[x] by 
collapsing the order parameter data for various different 
sizes, using the exponents j3/v = 1/5 and D = 5/2, which 
agrees perfectly with the MF prediction of the exponent 
values. Figure [3] is the scaling plot for the ER random 
networks with the noisy frequency distribution. As ex- 
pected, the scaling collapse is found with the same values 
of the exponents /3/P = 1/5 and v = 5/2. Moreover, the 
scaling function f[x] is identical for both the noiseless 
and noisy distributions with the same critical coupling 
strength K c . Differences can be seen only in higher-order 
finite-size effects. 

We also report the scaling plot for the z-regular ran- 
dom network with the noiseless frequency distribution in 
Fig. 0J As the strength of neighboring correlations in 
quenched networks depends on the details of networks, 
the critical point in the z-regular network is located dif- 
ferently from that in the ER random network, which is 
estimated as K c = 0.365(5) from the Binder cumulant 
analysis. Nevertheless, the corresponding exponent val- 
ues of P/v and v are identical to those in the ER random 
network as shown in Fig. 21 which are against the MF 
prediction discussed in the previous section. This im- 
plies that the random link-disorder fluctuations in the 
quenched z-regular network generate sample-to-sample 
Gaussian fluctuations of the same kind in the ER ran- 
dom networks as well as in the globally coupled network 
with frequency-disorder fluctuations. However, the scal- 
ing function f[x] by itself is distinct from that in the ER 
networks. 
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FIG. 2: (Color online) Scaling plot of A for various sys- 
tem size A^ in the ER random networks with the noiseless 
frequency distribution. 
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FIG. 3: (Color online) Scaling plot of A for various system 
size N in the ER random networks with the noisy frequency 
distribution. 
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V. SUMMARY 

In this work, we considered a finite population of cou- 
pled random frequency oscillators on the ER and the z- 



FIG 4: (Color online) Scaling plot of A for various system 
size N in the z-regular random networks with the noiseless 
frequency distribution. 



regular random network with link-disorder fluctuations. 
In the globally coupled network without any link disor- 
der, the finite-size fluctuations arc determined crucially 
by the existence of disorder in the intrinsic frequency 
distribution. Under the presence of link disorder in ran- 
dom networks, we find that the finite-size fluctuations 
are universal with the anomalous FSS exponent v = 5/2, 
independent of the details of the degree distribution (if 
not strongly heterogeneous) and also independent of the 
frequency disorder. We note, however, that the scaling 
function is not universal, depending on the details of the 
network connectivity. 
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